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Abstract 

A method is suggested to obtain solutions of the various quantum op- 
tical Hamiltonians in the framework of the asymptotic iteration method. 
We extend the notion of asymptotic iteration method to solve the 2x2 
matrix Hamiltonians. On a particular case, eigenvalues of the Rabi and 
Rashba Hamiltonians are computed. The method presented here repro- 
duces a number of earlier results in a natural way as well as leads to a 
novel findings. Possible generalizations of the method are also suggested. 
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1 Introduction 



The solutions of the quantum optical Hamiltonians are an important theme in 
the existing literature [1, 2, 3, 4, 5]. There has been a great deal of interest in 
quantum optical models which reveal new physical phenomena described by the 
2x2 matrix Hamiltonians. These Hamiltonians appear at the different fields of 
the physics [6, 7, 8, 9, 10, 11]. It is well known that the interaction of a two-level 
systems with a radiation field has a long varied history and these problems are 
often modelled by using 2x2 matrix Hamiltonians. As an example a physical sys- 
tem which describes the optical and electrical properties of confined electrons in 
semiconductor quantum wells, quantum dots and quantum wires depends on the 
Rashba spin orbit coupling and its equation is given by 2 x 2 matrix Hamiltonian. 
The Rashba splitting has been observed in many experiments and it constitutes 
the basis of the proposed electronic nano structures [12]. Due to the practical 
and technological importance of these models, it is not surprising that various 
aspects have been studied both analytically and numerically. Such systems have 
often been analyzed by using numerical methods because the implementation of 
the analytical techniques does not yield simple analytical expressions. Remark- 
ably, exact solutions have not been thus far presented except for special cases even 
though it has been suggested that the problem may be solved exactly and their an- 
alytical treatments require tedious calculations [13, 14, 15, 16, 17, 18, 19, 20, 21]. 

In recent years much attention has been focused on asymptotic iteration 
method (AIM) [22, 23, 24, 25, 26, 27, 28]. This method reproduces exact so- 
lutions to many exactly solvable differential equations and these equations can 
be related to the Schrodinger equation. It also gives accurate results for the 
non solvable Schrodinger equation including sextic oscillator, cubic oscillator, 
deformed Coulomb potential etc., which are important in applications to many 
problems in physics. Encouraged by its satisfactory performance through com- 
parisons with the other methods, we feel tempted to develop AIM to solve matrix 
differential equations. In contrast to the solution of the Schrodinger equation by 
using AIM including Coulomb, Morse, harmonic oscillator, etc. type potentials, 
the study of the quantum optical Hamiltonians has not attracted much attention 
in the literature. Such Hamiltonians have been found to be useful in the study 
of electronic properties of semiconductors, quantum dots and quantum wells. 

The aim of this paper is to develop AIM for solving 2x2 matrix Hamiltonians 
and discuss their applications. Results of our procedure include the solutions of 
the Rabi [11] and Rashba [6] Hamiltonians. We provide a first step toward the 
extension of the technique to the solution of the various matrix Hamiltonians, 
whose spectrum can not be obtained exactly. 

The paper is organized as follows. In section 2, we develop the AIM to obtain 



2 



eigenvalues and eigenfunctions of the wide range of the matrix Hamiltonians. 
Section 3 is devoted to solve the Rabi Hamiltonian in the framework of the AIM. 
In section 4 we present the solution of the Rashba Hamiltonian. In section 3 and 4 
we also discuss the bosonisation of the physical Hamiltonians whose original forms 
are given as differential operators. We also present a procedure to transform the 
bosons in to Bargman-Fock space which is necessary to obtain first-order matrix 
differential Hamiltonians. The paper ends with a brief conclusion. 

2 Formalism of the asymptotic iteration method 
for matrix Hamiltonians 

The AIM is proposed to solve the second-order differential equations and the 
details can be found in [22] . In this section we systematically extend the method 
for the 2x2 first-order matrix differential equations. We begin by rewriting a 
first-order differential equation in the following matrix form: 

I<t>' = «O0 (1) 

where <p = [0i, 4>2\' ', two-component column vector uq is 2 x 2 matrix function 
and / is 2 x 2 unit matrix. Note that <fi is the function of x and <f)' is the first 
derivative with respect to x. Now, in order to obtain a general solution to this 
equation in the framework of the AIM we use the similar arguments given in [22]. 
More explicitly, the differential equation (1) can be written as the two coupled 
equations 

0i = ^001 + ^002; 02 = C O02 + ^001 (2) 

where do, bo, cq and do are elements of the matrix uq. It is easy to show that n th 
derivative of the <p\ and 02 can be expressed as 



0" = Oi0i + 6i0 2 ; 02 = Cl0 2 + dl01 

0i" = «201 + &202 5 <t>2 — C 202 + ^201 

0i n) = O n _i0i + 6„_10 2 ; 02^ = C n _i0 2 + dn-101 

0i™ +1) = a«0i + bnfo; 02 n+1) = c n 2 + d n <f>i. 



(3) 



In order to discuss the asymptotic properties of (1), it is necessary to determine 
the coefficients a n , b n , c n and d n . After some straightforward calculation one can 
obtain the following relations: 

a n = aoa n _i + a' n _i + d b n _i 
b n = b Q a n ^i + b' n ^ + c 6„_i 

= c c n _i + d n _ x + Mn-l (4) 



On 



d n — d c n -i + d' n _ l + aoC?„_i. 
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Our task is now to introduce the asymptotic aspect of the method. For this 
purpose n th and (n + l) th derivative of the 0i and 2 can be written as 

A n ) (a i b n -i 
0i = a n ~i I 0i H 

0S n+1) = a„(0i + ^0 2 j; 0r l> = c n \^ 2 + ^0i J . (5) 

The coefficients a , b , d and c include the coupling constants. Therefore, for 
sufficiently large n we can suggest the following asymptotic constraints: 

71; = — = 72- (6) 




Q"n— 1 C n _i 

In this formalism the relations given in (6) imply that the wave function 0i and 
02 are truncated for sufficiently large n and the roots of the relations (6) belong 
to the spectrum of the matrix Hamiltonian. Therefore one can easily compute 
the eigenenergies of the Hamiltonian by solving (6) for the energy term when 
x — > as- 
under the asymptotic condition given in (6), one can find the wave functions 

0i and 02. When we take 1 (n) and 2 (n) by using (5) under the constraints 

given in (6), we obtain: 



0j n " ) = exp ( / dx ] or 2 n ' ) = exp ( / -^—dx ] . (7) 
\J a n _i J ' \J c n -! J 

Substituting the expression of a n (and c„) given in (4) into (7) and then replacing 
the 0i™' ) (and (p^) in the first expression of (5), respectively, one gets the following 
expressions; 

0i + 7i02 = exp (^J (a + 7id ) dx^j or 2 + 7 2 0i = exp (^J (c + 72&0) dxj . 

(8) 

An immediate practical consequence of these results is that the eigenvalues 
and eigenfunctions of the various quantum optical Hamiltonians can easily be 
determined. We mention here that there exist numerous physical Hamiltonians 
which can be written in the form of the first-order matrix differential equation 
or can be transformed in the form of the first-order matrix differential equation. 
For example, the Dirac equation [25] is a first-order matrix differential equation 
while most of the Hamiltonians of the quantum optical systems can be written as 
first-order matrix differential equation in the Bargman-Fock space. Now we can 
determine the eigenvalues and eigenfunctions of the 2x2 matrix Hamiltonians. 
In the following sections, it will be shown that this asymptotic approach opens 
the way to treatment of a large class of matrix Hamiltonians of practical interest. 
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3 Solution of the Rabi Hamiltonian by using 
AIM 



In this section we take a new look at the solution of Rabi Hamiltonian through 
the AIM. The Rabi Hamiltonian is a successful model of the interaction between 
matter and electromagnetic radiation. The Hamiltonian describes interacting a 
dipole interaction with a single mode of radiation in a two-level system. This 
Hamiltonian can be solved exactly within the framework of rotating wave ap- 
proximation and its quasi-exact or isolated exact solutions can be obtained for 
some specific cases. Let us consider the Rabi Hamiltonian for a two-level atom 

[11]: 

H — — — I + -mu 2 x 2 I + u cr + k (a + + a) x (9) 

2m ax 1 2 

where o;o is the atomic level splitting, uj is the frequency of the oscillation and k 
is the coupling strength of the atom to the field and, a , a + , ct_ are 



a 0=\ n 1 = n n h a - = 1 n ( 10 ) 






Pauli-related matrices. Eigenvalues of the Hamiltonian (6) can be obtained from 
the eigenvalue equation H<p = E(p where is a two component wave function. 
In general, it is difficult to determine the asymptotic aspects of the (9) in the 
present form. Therefore it is worth to transform (9) to an appropriate form. One 
way to transform a Hamiltonian in the form of the first-order matrix differential 
equation is to construct its bosonic representation. We are interested in the 
two-level system in one and two-dimensional geometry, whose Hamiltonians are 
given in terms of bosons-fermions or matrix-differential equations. Therefore, it 
is worth to express a suitable differential realizations of the bosons. Consider the 
following boson realizations 




a=\— — f- + J^-x; a + = -\l— — f- + J'-^x. (11) 
V 2mudx V 2// ' ! " ' ; ' 

They satisfy the usual commutation relation 

a,a + ] = 1; [a, a] = [a + ,a + ] = (12) 
In terms of the boson operators, the Hamiltonian (9) can be written as 



H = hu (a + a + ^] I + cu a + \\-^—k (o+ + cr_) (a + + a). (13) 
V 2/ V 2mcj 

It is interesting that this type Hamiltonian gives good results when it is solved by 
using the AIM in the Bargman-Fock space. In order to transform the Hamiltonian 
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to the Bargman-Fock space, we introduce the following transformation operator: 



r 1 = exp[^(a 2 + a +2 )]. (14) 

The action of the operator on the bosons is given by 

Tiar^ 1 = a cos ol\ — a + sin ol\\ TiG^r^ 1 = a + cos a\ + a sin a\. (15) 
For the values at\ = 7r/4 boson operators take the form 



1/ , \ h d , 1/ ,\ Irmo 

a ^T2 (a ~ a ) = V^ ; a ^71 (a + a WlT x (16) 

The final form of the Hamiltonian (13) with the realization (16) can be expressed 
as 

E = ^{*i + i) I+ ™ + fr°+ + '-){* + iL£)- (17) 

It is obvious that the Hamiltonian (17) can easily be written as the two coupled 
first-order differential equation 

(f)[ = a o 0i + b (f) 2 

02 = C O02 + d (j)i. (18) 

For the sake of simplicity let us take ft = w = m = 1 at that point, then the 
coefficients of the coupled differential equations are given by 



a 
b 



x (2E - 1 + k 2 - 2u; ) 

2x 2 - K 2 
k (1 - 2E - 2x 2 - 2u ) 

V2(2x 2 - k 2 ) 
x (2E - 1 + k 2 + 2uj ) 



« = g ... (19) 



K 



do = 



k (1 - 2E - 2x 2 + 2oj ) 
V2 (2x 2 - k 2 ) ' 



Using a simple MATHEMATICA program one can compute a n , b n , c n and d n 
by the relations given in (4). On the other hand, for each iteration, the quantiza- 
tion condition 8\{x) = b n -i(x)a n (x) — a n _i(x)b n (x) (and ^(x) = d n -i(x)c n (x) — 
c n-i{x)d n {x)) depends on some variables namely E u ,k,uj and x. It is noticed 
that the iterations should be terminated by imposing the condition Si(x) — 0, i — 
1, 2 as an approximation to Eq. (6) to obtain the eigenenergies. The calculated 
eigenenergies E n by means of this condition should, however, be independent of 
the choice of x. The choice of x is observed to be critical only to the speed of 
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the convergence of the eigenenergies, as well as for the stability of the process. 
In our study it has been observed that the optimal choice for x is the extremum 
point of the potential that is when x — 0. Therefore, we set x = at the end 
of the iterations. Then, the roots of the iteration produce eigenenergies for the 
Rabi Hamiltonian. The results are reported in Table I. 

We can also obtain eigenfunctions of the Hamiltonian (9) by using the rela- 
tion (8). Analytical expression for the ground state wave function of the Rabi 
Hamiltonian for some different parameters are given by 

k = 0; Uo = 0; 

1 

K = 0; Uo = -; 
k = i; u = 0; 

The solution of this system describes a quantum mechanical state of H provided 
that belongs to the Bargman-Fock space. It is noticed that when k or uoq is 
zero then the Hamiltonian (9) can be solved exactly. We have observed that for 
n — 5, accurate eigenvalues can be obtained after 15 iteration. As a consequence 
we have demonstrated that the solution of the Rabi Hamiltonian can be treated 
within the AIM. Our approach is relatively simple and gives accurate result. 



E = n + -; 



E = n; 



E = n + -; 



~~ \ C 2 x n ) ' 
C 2 x n ) ' 

_ / Cie~^(l + 2 v / 2x) n 
^ C 1 e~^(l + 2V2x) n 



(20) 



4 Solution of the Rashba Hamiltonian by using 
AIM 

The origin of the Rashba spin-orbit coupling in quantum dots is due to the lack of 
inversion symmetry which causes a local electric field perpendicular to the plane 
of heterostructure. The Hamiltonian representing the Rashba spin orbit coupling 
for an electron in a quantum dot can be expressed as [6] 

H R = ^ {Py(T x ~ P x (J y ) (21) 

where A# represents the strength of the spin orbit coupling, and it can be adjusted 
by changing the asymmetry of the quantum well via external electric field and 
the matrices a x , and a y are Pauli matrices. Now we assume that the electron is 
confined in a parabolic potential 

V = l -m^l{x 2 + y 2 ) (22) 

here m* is the effective mass of the electron and uo Q is the confining potential 
frequency. The Hamiltonian describing an electron in two-dimensional quantum 
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dot takes the form 



H = ( P * + p y) + ^Ba, + V + H R . (23) 

The term ^gfiBa introduces the Zeeman splitting between the (+)x— polarized 
spin up and (—)x— polarized spin down. The factors g is the gyromagnetic ratio 
and \x is the Bohr magneton. The kinetic momentum P = p + eA is expressed 
with canonical momentum p = —ih(d x , d y , 0) and the vector potential A can be 
related with the magnetic field B = V x A. The choice of symmetric gauge vector 
potential A = B/2{— y, x, 0), leads to the following Hamiltonian 

h 2 ( d 2 d 2 \ 1 „ 2 , 2 2 ^ 

+-ihu c [y— - x—J + -gfiBco + H R (24) 

where u c = eB/m* stands for the cyclotron frequency of the electron, uj = 

\Jujq + is the effective frequency. From now on we restrict ourselves to 

the solution of (24). 



4.1 Bosonisation of the Hamiltonian 

The Hamiltonian (24) can not be solved within the framework of the AIM in 
the present form. Our task is now to demonstrate that the Hamiltonian (24) 
can be expressed as two coupled first-order differential equation in the Bargman- 
Fock space. One way to express the Hamiltonian H with boson operators is 
to introduce an appropriate differential realization for bosons. It can easily be 
bosonised when the boson operators are realized as 



b+ = i/^-^-te^-^ (25) 



b = ^ — (x + iy) + ] l — (d x + id y ). 

They satisfy the usual commutation relations. Insertion of (25) into (24) yields 
the following Hamiltonian: 

H = ^(a + a + 6 + 6 + l) + ^(a+a-6+6) 
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h 



-A 



R 



(b + - a)a + + (b - a + )a_ 



(26) 



Now we turn our attention to the transformation of the Hamiltonian in the form 
of the first-order one variable matrix differential equation. A simple connec- 
tion between the Hilbert space and the Bargman-Fock space can be obtained by 
transforming the differential realizations of the creation and annihilation opera- 
tors (25). This can be done by introducing the following similarity transformation 
operators 



cxp 



^ (b 2 + 6 +2 )| ; A = exp [(3 (a + b + b + a) 



(27) 



The operators act on the bosons as follows: 



rvrr 1 = 

r 2 6r 2 - 1 = 

AaA- 1 = 

AbA- 1 = 



acos«i - a sm«!; F 1 a + F^ — a + cos ot\ + a sin ot\ 
a cos a 2 — a + sin a 2 ; I^c^r^ 1 = a + cos a 2 + a sin a 2 
a cos (3 — b sin f3; Aa + A _1 = a + cos (3 — b + sin (3 
b cos (3 + a sin (3; A6 + A _1 = b + cos (3 + a + sin (3 



(28) 



These transformations play a key role to construct one variable first-order matrix 
differential equation form of (26). The similar transformation and the change of 
the variable y — > iy gives the following realizations 



Ariar^A -1 = 



Arw^A- 1 



/ h d 

m*u dx ' 



/ h d 
m*uj dy : 



b + 



> AT 1 a+T^ 1 A- 1 = 
ATzb+T^A' 1 = 4 



irrruj 



-x 



h 



y- (29) 



It is obvious that the Hamiltonian (26) can be put in the form of the first- 
order matrix differential with the realization of bosons (29), but it consists of 
two variables x and y. In order to separate the variables consider the following 
conserved quantity of the Hamiltonian (26): 



K = a + a — b + b an. 

2 



(30) 



If K and H commute, the eigenfunction of K is also eigenfunction of the H. 
Therefore it is worth to obtain eigenfunction of K. When we solve the following 
eigenvalue equation 

K |ni n 2 ) = (k + jn x n 2 ) , (31) 
in the Bargman-Fock space we obtain the following expression: 



il>(x, y) = x k <j) (xy) | T) + x k+1 <j) (xy) 1 1) . 



(32) 
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where |f) stands for up state and ||) stands for down state. The eigenfunction 
of the Hamiltonian can be obtained from the relation 



\n!,n 2 ) = T^T^A 1 i>(x,y). 



(33) 



Substitution of (32) into the Hamiltonian (26) by using the realization of bosons 
in (29) leads to the following set of one variable coupled differential equations 



hu 



~ d , ku c 
2z— + k + l + —^ 

dz 2u 



+A 



R 



»gB E 
2hu hu 

m*uz d 
fc + 1 — + z— 

n az 



01 (Z) 

02 (Z) 



(34) 



and 



hu 



2z— + k + 2 + 

az 



2u 



+ -^z ir- 



+A 



R 





E 


2hu 


hu 


m*u 


d 


h 


dz 



02 (Z) 

01 (z) = 



(35) 



where z = xy and E is the eigenvalues of the Hamiltonian H and 0i (z) and 
02 (z) correspond up and down eigenstates of the Hamiltonian H, respectively. 
Following the analysis of [8] which was constructed to obtain the solution of the 
Rashba Hamiltonian, one can obtain the quasi-exact solution of the differential 
equations. Here we present solution of the problem in the framework of the AIM. 



4.2 Solution of the Rashba Hamiltonian 

In the previous section we have formulated the Rashba Hamiltonian based on the 
two boson operators and we have discussed its transformation to the one variable 
differential equation. In this section we present the AIM which leads to solution 
of the Hamiltonian (26). For the sake of simplicity let us take K — m* — u — 1, 
then the last expression given above can easily be written as the two coupled 
differential equation 

01 = OO01 + ^002 

02 = C O 02 + d (j)i (36) 

where 

u (2E + X 2 R + Bg/2 - u c k - 2u(k + 1)) 
a ° ~ (4u*z + Al) 

X r (-2E + Bgn + u c (k + 1) - 2uk + 4u 2 z) 
~ 2 {4u 2 z + A|) 
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uz (2E - Bg/2 - u c (k + 1) - 2uj(k + 2)) - A|(l + k - uz) 

Co " ; {a^z + as (37) 

A fi (2£ + - u c k - 2u(k + 1) - 4cj 2 z) 
~ 2z {4u 2 z + A|) 

As in the previous section, one can compute a n ,b n , c n and d n using the relations 
given in (4) by a simple MATHEMATICA program. When z — 0, the eigenvalues 
are produced by the roots of the quantization condition. The computed results 
of the iteration are reported in Table 2. 

As it is seen that the energy eigenvalue equation is easily obtained by using 
AIM. This is the advantage of the AIM that gives the eigenvalues directly by 
transforming the quantum optical Hamiltonians in the form of the 2x2 matrix 
differential equation form. In a similar manner as in the previous section one can 
also obtain eigenfunction of the Hamiltonian (34) and (35) by using the relation 
(8). 



5 Conclusion 

We study the AIM to solve the problem of an electron in a quantum dot in the 
presence of both magnetic field and spin-orbit coupling. Our formulation gives 
an accurate result for the eigenvalues of Rabi and Rashba Hamiltonians. The 
suggested approach can easily be modified to solve some other quantum optical 
problems. 

Furthermore we have presented a transformation procedure that offers several 
advantageous, especially if one wishes to describe the eigenvalues of the bosonic 
Hamiltonians by using AIM. It is obvious that the technique presented here have 
been used in a variety of problems to compute their spectrums. We have presented 
the steps towards an extension of the AIM. 

The technique given in this article can be extended in several ways. The 
Hamiltonian of a quantum dot including position dependent effective mass may 
be formulated and solved within the procedure given here. We hope that our 
method leads to interesting results on the spin-orbit effects in quantum dots in 
future research. Success of our analysis leads to the solution of the wide range of 
physical systems. 

One of the authors (O. Ozer) is grateful to the Abdus Salam International 
Centre for Theoretical Physics, Trieste, Italy, for its hospitality. 
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Table 1: Eigenvalues of the Rabi Hamiltonian for different values of energy state 
n, splitting constant uq, and coupling constant k. 



n 


u 


K=0.0 


/t=0.25 


k=0.50 


k=0.75 


«=1.0 




0.0 


0.5 


0.46875 


0.375 


0.21875 


0.0 






0.5 


0.46875 


0.375 


0.21875 


0.0 





0.5 


0.0 


-0.015748 


-0.064513 


-0.151195 


-0.284922 






1.0 


0.8080812 


0.5865567 


0.3344402 


0.0488330 




1.0 


-0.5 


-0.510489 


-0.542870 


-0.600162 


-0.688478 






0.5 


0.4494499 


0.3117158 


0.1101679 


-0.141057 




0.0 


1.5 


1.46875 


1.375 


1.21875 


1.0 






1.5 


1.46875 


1.375 


1.21875 


1.0 


1 


0.5 


2.0 


1.7359493 


1.4482320 


1.1416351 


0.8206445 






1.0 


1.1601606 


1.2791925 


1.3249343 


0.0488339 




1.0 


1.5 


1.4124109 


1.1964923 


0.9129976 


0.5913838 






1.5 


1.5294031 


1.5980360 


1.6078394 


1.2925094 




0.0 


2.5 


2.46875 


2.375 


2.21875 


2.0 






2.5 


2.46875 


2.375 


2.21875 


2.0 


2 


0.5 


2.0 


2.2317488 


2.4052487 


2.4171516 


2.1405752 






3.0 


2.6813089 


2.3498342 


2.0376015 


1.8229775 




1.0 


2.5 


2.5660788 


2.6980532 


2.5296180 


2.0954608 






2.5 


2.3778828 


2.1020231 


1.8333934 


1.8464607 




0.0 


3.5 


3.46875 


3.375 


3.21875 


3.0 






3.5 


3.46875 


3.375 


3.21875 


3.0 


3 


0.5 


4.0 


3.6359015 


3.2763600 


3.0202409 


2.9667546 






3.0 


3.2858045 


3.4872316 


3.4040837 


3.0257950 




1.0 


3.5 


3.3455073 


3.0291733 


2.9681181 


3.1101317 






3.5 


3.6001795 


3.7530103 


3.4089357 


2.9237568 




0.0 


4.5 


4.46875 


4.375 


4.21875 


4.0 






4.5 


4.46875 


4.375 


4.21875 


4.0 


4 


0.5 


4.0 


4.3305823 


4.5395425 


4.347928 


3.9264039 






5.0 


4.5965321 


4.2235326 


4.070904 


4.0877781 




1.0 


4.5 


4.6320470 


4.7466147 


4.294056 


3.7849395 






4.5 


4.3150285 


3.9930696 


4.121917 


4.2669272 




0.0 


5.5 


5.46875 


5.375 


5.21875 


5.0 






5.5 


5.46875 


5.375 


5.21875 


5.0 


5 


0.5 


6.0 


5.5615693 


5.1903384 


5.1469657 


5.1455854 






5.0 


5.3692738 


5.5672925 


5.2808323 


4.8635200 




1.0 


5.5 


5.2862655 


5.0134462 


5.2558821 


5.2799593 






5.5 


5.6619147 


5.7022888 


5.1889935 


4.7100710 
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Table 2: Eigenvalues of the Rashba Hamiltonian for different values of external 
magnetic field magnitude, B, energy state n, and coupling constant k. 



k 


B 


n 


K = 0.0 


k = 0.25 


k = 0.5 


k = 0.75 


K = 1.0 









0.0 


0.0 


0.0 


0.0 


0.0 






1 


1.0 


0.9393054 


0.7737872 


0.5320678 


0.2330309 







2 


2.0 


1.9962478 


1.9465668 


1.7810714 


1.4877948 






3 


3.0 


2.9429223 


2.8214396 


2.7072634 


2.5694541 






4 


4.0 


3.9928041 


3.9101840 


3.6989847 


3.4345727 







5 


5.0 


4.9462371 


4.8526593 


4.7537468 


4.5359330 









1.5 


1.5 


1.5 


1.5 


1.5 






1 


0.5 


0.4776216 


0.4097986 


0.2944479 


0.1280664 




3/2 


2 


3.0 


2.8877702 


2.6235301 


2.2882565 


1.9132967 






3 


4.0 


4.0892712 


4.2713119 


4.2790928 


3.8467145 






4 


5.5 


5.3107639 


4.9316662 


4.6839650 


4.8228044 






5 


6.5 


6.6655604 


6.9171192 


6.5731725 


6.0024178 









-1.0 


-1.0 


-1.0 


-1.0 


-1.0 






1 


2.0 


1.8847720 


1.6082091 


1.2505849 


0.8439143 







2 


3.0 


3.0463969 


3.0303266 


2.7527743 


2.3150943 






3 


4.0 


3.8955162 


3.7596438 


3.7797919 


3.6831314 






4 


5.0 


5.0363268 


4.9398058 


4.6186235 


4.3781510 






5 


6.0 


5.9052818 


5.8399949 


5.8328977 


5.5201455 


1 







1.0 


1.0 


1.0 


1.0 


1.0 






1 


2.5 


2.4555260 


2.3240428 


2.1107652 


1.8227033 




3/2 


2 


5.0 


4.8690809 


4.5668665 


4.1837412 


3.7484901 






3 


6.0 


6.1078281 


6.3230999 


6.2570599 


5.7506451 






4 


7.5 


7.2943805 


6.8910658 


6.7283380 


6.9481530 






5 


8.5 


8.6817766 


8.9384866 


8.5213644 


7.9201055 
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